Topologically protected flat zero-energy surface bands in non-centrosymmetric superconductors 
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Nodal non-centrosymmetric superconductors (NCS) have recently been shown to be topologically non-trivial. 
An important consequence is the existence of topologically protected flat zero-energy surface bands, which are 
related to the topological characteristics of the line nodes of the bulk gap via a bulk-boundary correspondence 
[Schnyder and Ryu, arXiv: 101 1.1438 1. In this paper we examine these zero-energy surface bands using a 
quasiclassical theory. We determine their spectrum and derive a general condition for their existence in terms 
of the sign change of the gap functions. A key experimental signature of the zero-energy surface bands is a 
zero-bias peak in the tunneling conductance, which depends strongly on the surface orientation. This can be 
used as a fingerprint of a topologically non-trivial NCS. 
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Introduction. A key experimental signature of topologi- 
cal insulators and superconductors is the existence of topo- 
logically protected zero-energy surface or edge states, some 
of which are of Majorana typel^But zero-energy boundary 
modes of topological origin can also occur in gapless topolog- 
ical systems that exhibit topologically stable Fermi points or 
in nodal superconductors with non-trivial topology. In such 
systems one generically finds topologically protected disper- 
sionless zero-energy surface states, i.e., flat bands at the sur- 
face. Such flat bands are known to occur at the zig-zag and 
bearded edge in graphene, 4 5 on the (110) surface of a d x 2_ y 2- 
wave superconductor,^ within vortices in the A phase of 
3 He,^and in other systems with topologically protected Dirac 
points P 

It has recently been realized that nodal non-centrosymme- 
tric superconductors (NCS) are topologically non-trivial states 
of matter. 10 12 An NCS is realized in a system lacking inver- 
sion symmetry, which gives rise to an antisymmetric spin- 
orbit (SO) coupling and consequently to the admixture of 
spin- singlet and spin-triplet components in the superconduct- 
ing state. There is a steadily growing list of these remark- 
able materials, most notably Li 2 Pd^Pt 3 _ x B!EEl YaCsP^and 
the heavy-fermion compounds CePtsSip^ CeRhSispl and 
CelrSis It is predicted that the topological non-triviality of 
a nodal NCS leads to topologically protected zero-energy sur- 
face bands, which only occur within regions of the surface 
Brillouin zone bounded by the projected nodal lines of the 
bulk gap. Since these zero-energy surface bands give a singu- 
lar contribution to the surface density of states, we can expect 
them to lead to a zero-bias conductance peak (ZBCP). 

It is the aim of this paper to investigate the appearance of 
zero-energy surface bands in an NCS using the quasiclassi- 
cal scattering theory. 7 This method is ideal for exploring the 
bound surface states of unconventional superconductors, and 
has revealed key aspects of the surface physics of the cuprate 
high-T c compounds 67 and NCS systems. 19 24 We hence de- 
rive the surface-bound- state spectrum and a general condition 
for the existence of the zero-energy surface bands in terms of 
a sign change of the superconducting gap function across the 
Fermi surface. This condition is complementary to the topo- 
logical criterion given in Ref . Q21 We then compute the tun- 



neling conductance between a normal metal and an NCS as a 
function of both the surface orientation and the relative mag- 
nitude of the spin-singlet and spin-triplet pairing states. We 
argue that the strong dependence of the ZBCP on both these 
variables provides a powerful diagnostic test of the pairing 
state of an NCS. 

Model system. — We phenomenologically model the NCS as 
a single-band system, described by the Bogoliubov-de Gennes 
(BdG) Hamiltonian 

«(W-f e(k) + g(k) -* A(k) ^ m 
n W-{ At(k) -,(k)+g(k).<7*j' (1) 

where e(k) is the spin-independent part of the band disper- 
sion, g(k) = — g(— k) is the antisymmetric SO coupling, and 
A(k) = (A s + d(k) • a) (icr y ) is the superconducting gap 
function. It is convenient to express the Hamiltonian Eq. ([T]) 
in the so-called helicity basis which diagonalizes the kinetic 
term, yielding two helicity bands with dispersions ^(k) = 
£(k) ± |g(k)|. In the absence of interband pairing the criti- 
cal temperature is maximized by taking the spin-triplet pairing 
vector d(k) to be aligned with the polarization vector of the 
SO coupling g(k), 25 i.e., we parametrize the triplet compo- 
nent of the gap function and the SO coupling as d(k) = A £ l k 
and g(k) = cd k , respectively. Hence, the gaps on the two he- 
licity bands are = A s ± A t |l k | = A (q ± |l k |)/(<? + 1), 
where the parameter q = A s /A t interpolates between pure 
triplet pairing (q = 0) and pure singlet pairing (q — >> oo). For 
simplicity we assume the pairing amplitudes A s and A t to be 
constant and have positive sign. 26 We note that higher-order 
angula r momentu m components of the gap A have also been 
studied! 21 ' 23 ' 29 ' 30 ' 

The specific form of the pseudovector l k depends on the 
symmetries of the non-centrosymmetric crystal P 7 Ignoring 
the periodic Brillouin-zone structure we employ a small- 
momentum expansion which for the tetragonal point group 
C± v (relevant for CePtaSi, CeRhSi3, and CelrSis) gives, to 
lowest order, the symmetry- allowed form 

l k = k y ± - k x y. (2a) 

For the cubic point group O, which is relevant for 
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FIG. 1. (color online) Panels (a), (c) and (e): Variation of the sign of 
A^ over the Fermi surface as seen along the [111] direction for (a) 
the Cav point group with q = 0.5, (c) the O point group with q — 
0.35, and (e) the Td point group with q = 0.4. Dark red indicates 
sgn(A^) = 1, whereas light blue is sgn(A^) = —1. Panels (b), 
(d) and (f): Surface bound states at the (111) face corresponding to 
the same parameters as in (a), (c) and (e), respectively. The color 
scale of each plot is normalized by the maximum bound state energy, 
(b) E m ax = 0.333 A , (d) E ma x 

= 0.211 A , and (f) 
0.421 Ao. White space indicates the absence of a bound state, and 
the circle is the projection of the Fermi surface, i.e., |k|| | = Jzf- 



^PcLpPta-^B, the expansion of lk reads 

Ik = k x (l + g 2 [k 2 y + k 2 z }) x + k y (1 + g 2 [k 2 x + k 2 ]) y 

+k z (l + g 2 [kl + kl])z. (2b) 

For the spherical Fermi surfaces considered below, the gap 
A^ only has nodes if g 2 ^ 0. Finally, we also consider 
the tetrahedral point group Td, experimentally represented by 
Y2C3, where lk takes the form 

lk = k x (k 2 y - k 2 z )± + k y (k 2 z - kl)y + k z (k 2 x - k 2 y )z. (2c) 

Surface bound states. — We wish to solve the BdG equa- 
tion H^(r) = E^(r) at a surface for states which decay into 



the bulk, i.e., which have energies lying within the gap. We 
define coordinates parallel (ry) and perpendicular (r±) to the 
surface such that the NCS occupies the half-space r± > 0. 
For simplicity, we ignore the SO splitting of the bands and as- 
su me coin cident spherical helical Fermi surfaces with radius 
fcp ! 2Q l 21 l 24 l We have verified that relaxing this approximation 
does not alter the condition for the existence of zero-energy 
surface bands. It is therefore convenient to regard the momen- 
tum components in the definition of lk [Eq. ([2])] as normalized 
by kp> If we take g 2 = —1.5 in Eq. (2b), then for all point 
groups we find that A^ has point nodes for q = 0, 1, line 
nodes for < q < 1, and is fully gapped for q > 1. Examples 
of the nodal structure of A^ are shown in Fig. [I] 

We obtain the following ansatz for the surface bound state 
wavefunction 



(3) 



: p k.k 



where k = (k|i,fej_) and k = (kn,— k±) are wavevectors 
with |k| = |k| = kp- The momentum component parallel 
to the surface, ky, is a good quantum number due to transla- 
tional invariance. The positive and negative helicity compo- 
nents ^±(p) are given by 
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respectively, with 7p = (Ap) 1 [E — isgn(p±)(\A^\ 2 — 
E 2 ) 1 / 2 ]. The wavefunction components decay into the bulk 
over the inverse length scale = m(h 2 \kj_\)~ 1 (\A^\ 2 — 
E 2 ) 1 / 2 , where m is the effective mass. 

A bound state occurs if the coefficients a n (p) in Eq. ^ 
can be chosen so that the wavefunction vanishes at the sur- 
face. After some algebra, this yields the following condition 
for surface bound state formation, 



= ( 7 ±- 



7k)(7 £ -7 k + )(|lk||lk|- 
"7£)(7~ -7k)(|lk||lkl 
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(5) 



Setting E = in Eq. ^ and observing that 7p | £ „ 
— zsgn(px)sgn(Ap ), we find that Eq. ^ has a non-trivial 
zero-energy solution whenever (i) sgn(A^) = sgn(A~ ) = 



-1 and l k • lr 



|l k ||l £ | or (ii) sgn(A-) 



-sgn(Ar). 



The latter condition corresponds to the topologically protected 
zero-energy surface bands found in Ref. [T21 which occur 
within a finite region of the surface Brillouin zone bounded 
by the projected nodes of the bulk gap. It is interesting that 
this condition is quite similar to the one for zero-energy states 
at surfaces of unconventional centro symmetric superconduc- 
tors, i.e., the sign of the gap must reverse between k and k. It 
is therefore somewhat surprising that the relative sign between 
the gaps on different helicity bands at these wavevectors is ir- 
relevant, despite the mixing of the helicity components in the 
edge- state wavefunction. 



In Figs. [TJb), (d), and (f) we plot the dispersion of the 
bound states appearing at the (111) surface of an NCS with 
point group C^, O, and Td, respectively. For all three point 
groups we find zero-energy surface bands for < q < 1, 
i.e., there are two-dimensional regions in the surface Brillouin 
zone where the energy of the bound states vanishes [black re- 
gions in Figs.[lTb), (d), and (f)]. As can be seen by compar- 
ison with Figs7[lja), (c), and (e), these occur where the sign 
of on the forward-facing half of the Fermi surface is dif- 
ferent to that on the backward-facing half. This is a manifes- 
tation of the bulk-boundary correspondence found in Ref . [12] 
which relates the topologically stable nodal lines of the bulk 
gap to the zero-energy surface bands. A quantized topologi- 
cal invariant is associated with the line nodes, which protects 
them and hence also the zero-energy surface bands from suf- 
ficiently small symmetry preserving perturbations. Large per- 
turbations which remove the nodal rings (e.g., increasing the 
singlet-to-triplet ratio q past 1) therefore also destroy the as- 
sociated zero-energy surface bands. 

The zero-energy surface bands are in general accompanied 
by dispersing modes [colored regions in Figs.[T](b), (d), and 
(f)], which occur in regions of the surface Brillouin zone 
where sgn(A^) = sgn(A~) = —1. The coexistence of 
dispersive and non-dispersive surface states give rise to intri- 
cate bound-state spectra. In the case of the C^ v point group 
[Fig. [ija)] the linearly dispersing states at sufficiently small 
| &2, || I ar i se f rom me same mechanism as the states found along 
the (100) direction in Ref. [20l Note the line of zero-energy 
states at k\ n =0 connecting the two flat bands. 

Tunneling conductance. — Although the surface bound 
states do not form at an interface with a normal metal, for 
a low-transparency barrier between a metal and an NCS, the 
physical mechanism discussed above leads to the formation of 
interface-resonance states. For a surface which displays zero- 
energy bands, it is well-known that the corresponding tunnel 
junction will show a sharp ZBCP in the low- temperature tun- 
neling conductance. As in the case of the d-wave gap in the 
high-T c cupratesP^J a direction-dependent ZBCP in an NCS 
due to the topologically protected zero-energy surface bands 
would be a key experimental signature of the pairing state. 

The zero-temperature charge conductance as(eV) for tun- 
neling into the NCS is a ge neralizati on of the usual Blonder- 
Tinkham-Klapwij k 

formula 19 2 2EH2l 



as (eV) = J2U+\i: 

ky [ <J,(J> 



, (6) 



where aZ' a and bZ. ,CT are the Andreev and normal reflection 
for electron injection into the NCS, respectively. These coef- 
ficients are determined by solving the BdG equation for the 
junction at energy E = eV, with the wavefunction ansatz 
Mk| h r) = ©(-^)^k„,a( r ) + 0(^k |hCT (r) where 
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FIG. 2. (color online) Panels (a), (c) and (e): Normalized conduc- 
tance spectra at the (111) face for an NCS with point group (a) Ca v , 

(c) O, and (r) Td and various values of < q = A s / A t < 1. In all 
panels we take Z = 3 and T = K. Curves are vertically shifted by 
multiples of 0.2. Panels (b), (d) and (f): ZBCP as a function of sur- 
face orientation for an NCS with point group (b) Cav with q = 0.5, 

(d) O with q = 0.35, and (f) T d with q = 0.4. The color at each 
point on the sphere indicates the height of the ZBCP for the corre- 
sponding normal vector. Panel (g): Variation of ZBCP height at the 
(111) face as a function of q. 
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■dj;>n(k)e ik - r ] . (7) 



^e,a = U 1 + M - ^°>°) T and t/> hi<T = 1(0,0,1 + 
<j, 1 — cr) T are the electron and hole spinors in the normal 
metal, respectively. We adopt the assumption that the bias 
energy is small compared to the Fermi energy so that the 
wavevectors in Eq. ([7]) are well- approximated to have mag- 
nitude fcpP^EH xhe insulating barrier at r±_ = is mod- 
eled as a c) -function of height U. The coefficients in Eq. ([7} 
are then chosen such that ^ a is continuous at the interface, 
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i.e., ^ a (k||,r)| r ^ =0 - = ^(ky, r)| r ^ =0 +, while the deriva- 
tive obeys d r± ^a(^\\, r)| r ^ =0 + - d r± ^ a (k^r)\ r±=0 - = 
2Z\I/ cr (k||, r)\ r±= o, where Z = mU/H 2 with m is the effec- 
tive mass, assumed the same in the NCS and the metal. 

In Fig. |2ja), (c), and (e) we show the conductance spectra 
at the (111) face for the C± v , O, and Td point groups, respec- 
tively. The spectra are normalized by the normal- state con- 
ductance (j 7v = Ek l k ll ?/( z<2 + l k || I 2 )- For < g < 1 we 
find a sharp ZBCP, signaling the existence of the zero-energy 
surface bands. For the most part, this peak is well-separated 
from the contributions of dispersing states or the edges of the 
bulk gap. A notable exception is the C± v case [Fig. [2ja)], 
where for small q we find the ZBCP superimposed on a wide 
dome-like feature due to the dispersing surface modes ES2I As 
shown in Fig.[2jg), the ZBCP vanishes as q — >> 1 and is absent 
for q > 1, consistent with the topologically trivial gapped 
state which forms when the singlet pairing dominates. We 
note that in the case of the O and point groups, the ZBCP 
height shows kinks as a function of q, marked by the arrows 
in Fig.|2jg). At these values of q there is a Lifshitz-type transi- 
tion in the BdG spectrum at which the nodal rings touch each 
other and then reconnect. 

Finally, we consider the direction dependence of the ZBCP 
[Fig.[2jb), (d), and (f)]. The strong variation in the conduc- 
tance with the surface orientation reflects the changing pro- 
jection of the bulk nodal lines onto the surface Brillouin zone. 
Note that for all point groups the ZBCP is absent along the 
crystal axes. This can be exploited as a test of the pairing state 
in an NCS: For example, the observation of a ZBCP along 



the (111) direction, but its absence along the (100) direction, 
would lend strong support to a model of an NCS pairing state 
with a dominant triplet component. 

Conclusions. — In this paper we have used quasiclassi- 
cal scattering theory to study the appearance of topolog- 
ically protected zero-energy bands at the surface of non- 
centro symmetric superconductors for three experimentally 
relevant choices of the point group. We have derived a general 
condition for the existence of these states, which is consistent 
with the bulk-boundary correspondence found in Ref.[l2l The 
surface-bound- state spectrum has been computed and shown 
to allow the coexistence of zero-energy flat bands with dis- 
persing states. We have also calculated the tunneling conduc- 
tance, where the presence of these zero-energy states mani- 
fests itself as a ZBCP. The ZBCP displays a strong depen- 
dence on the surface orientation and in particular vanishes 
along the crystal axes. We propose that this dependance can 
be exploited as a test of the orbital and spin pairing symme- 
tries in these materials. 

Prospects for further work are promising. For example, 
the large degeneracy of the zero-energy surface bands may 
be expected to lead to instabilities towards symmetry-broken 
states in the presence of interactions. The possible coexis- 
tence of time-reversal- symmetry -breaking and time-reversal- 
symmetry-preserving order parameters near the surface of an 
NCS is particularly tantalizing. 
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